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Ab initio methods for electronic structure of molecules have reached a satisfactory accuracy for 
calculation of static properties, but remain too expensive for quantum dynamical calculations. We 
propose an efficient semiclassical method for evaluating the accuracy of a lower level quantum 
dynamics, as compared to a higher level quantum dynamics, without having to perform any quantum 
dynamics. The method is based on dephasing representation of quantum fidelity and its feasibility is 
demonstrated on the photodissociation dynamics of CO2. We suggest how to implement the method 
in existing molecular dynamics codes and describe a simple test of its applicability. 



Ab initio methods for electronic structure of molecules 
have reached a satisfactory accuracy for calculation of 
static properties, such as energy barriers or force con- 
stants at local minima of the potential energy surface 
(PES). The most accurate of such methods remain out 
of reach when one wants to describe molecular proper- 
ties depending on the full quantum dynamics [T]. To 
make a calculation feasible, one has to approximate the 
dynamics of the system [2J [3] or the PES pQ , but both 
approaches can have nontrivial effects on the result [I]. 
We consider only the second approach, in which quan- 
tum dynamics is done exactly but on a PES obtained by 
a lower level electronic structure method that is less ac- 
curate but also less expensive. When such a calculation 
is finished, its accuracy is not known because of the for- 
bidding expense of the dynamics on the more accurate 
potential. In this Letter, we propose an efficient and ac- 
curate semiclassical (SC) method for evaluating the accu- 
racy of the lower level quantum dynamics without having 
to perform the higher level quantum dynamics. Since our 
method does not even require computing the lower level 
quantum dynamics, it can also be used to justify, in ad- 
vance, investing computational resources into the lower 
level calculations. 

For simplicity, we use the Born-Oppenheimer approx- 
imation and focus on the quantum dynamics of nuclei, 
although our method applies to any quantum dynam- 
ics and should be valid even when nonadiabatic effects 
are important. The time-dependent Schrodinger equa- 
tion is solved in two stages: First, the time-independent 
equation for electrons is solved with fixed nuclear con- 
figurations, and then nuclear motion is calculated on the 
resulting electronic PES. We consider three PESs: T^xact 
is the exact PES that describes our system. V acc is a very 
accurate high-level electronic structure PES (presumably 
"almost exact"), which is too expensive to be used for 
quantum dynamics. V appr is an approximate PES, ob- 
tained by a lower-level electronic structure method (or 
by an analytical fit of V acc ) and "cheap" enough to be 
used for quantum dynamics. 

Various quantities can describe different time- 



dependent features of a quantum system, but a single 
quantity that includes all information about the sys- 
tem is the time-dependent wave function, yb(t). One 
way to evaluate the accuracy of quantum dynamics 
on the approximate PES would therefore be to com- 
pute the quantum-mechanical (QM) overlap /qm(^) := 
( -! /'cxact(£)|' i / ; appr(i))) where the subscript oiip denotes the 
corresponding PES used for propagating the initial state 
■0(0). The quantity F(t) :— |/qm(^)| 2 is known as quan- 
tum fidelity or Loschmidt echo, and has been defined by 
Peres [5] to measure the sensitivity of quantum dynamics 
to perturbations. Much effort has been devoted to the 
study of temporal decay of fidelity and many universal 
regimes have been found [B]. In our setting, if F(t) w 1 
for all times up to t max , we can trust quantum dynamics 
on the approximate potential V appI and use the resulting 
V'appr(i) to compute all dynamical properties up to t max - 
In calculations, we do not know T4 xact and must use the 
accurate potential V acc , and so we approximate /qm by 

/QM := (0acc(*)|0appr(t)) • (1) 

Since V acc is too expensive, ipa, cc {t) cannot be computed. 
We describe a method which gives an accurate estimate 
of /qm without having to compute ip&cc(t) nor f/'appr^)- 
The method is based on the dephasing representation 
(DR) of quantum fidelity, a SC approximation proposed 
by one of us to evaluate fidelity in chaotic, integrable, 
and mixed systems even in nonuniversal regimes sensi- 
tive to the initial state and details of dynamics [7J |8J. 
Presently, we are interested in a specific type of "pertur- 
bation," namely the difference AV = Va Ppr — V acc between 
the approximate and accurate PESs. The DR of fidelity 
amplitude is an interference integral 

/drW := / dx° pwix ) exp [-iAS(x°, t)/h] , (2) 

AS(x°,t) = fdrAV^x )]. (3) 
Jo 

Here x denotes a point (q,p) in phase space, the super- 
script is the corresponding time, AS(x°,t) is the action 
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due to AV along the trajectory ql cc (x°) of V acc , and pw 
is the Wigner function of the initial state ip, 

pw(x) = h- d J dtr (q+ ^ (« - f) cx p ( i ^r) • 

In "dephasing representation," all of fidelity decay ap- 
pears to be due to interference and none due to decay of 
classical overlaps. 

Surprising accuracy of the DR was justified by the 
shadowing theorem [7], or, in the case of the initial state 
ip(0) supported by a Lagrangian manifold, by the struc- 
tural stability of manifolds [9]. Interestingly, validity of 
the DR goes much beyond validity of the SC approxi- 
mation for quantum dynamics on V acc or V appT . Quali- 
tatively, this is due to mitigating the "sign problem" of 
quantum or SC dynamics: large actions needed for dy- 
namics on V acc or Va Ppr are replaced by much smaller 
actions AS needed for fidelity calculation. Rapid oscil- 
lations in the SC expression for the dynamics on V acc or 
T4ppr are replaced by much slower oscillations in the DR. 
In a chaotic system where ~ 10 35 trajectories would be 
needed for computing ip acc (t) semiclassically, as few as 
1000 trajectories were sufficient to compute fidelity am- 
plitude |10j . Accuracy of DR was explored numerically in 
Refs. [7HHHUJ which suggest that Fdr starts to deviate 
from Fqm after the Heisenberg time tn — h/AE where 
AE is the mean level spacing. Errors of DR can be esti- 
mated analytically, suggesting that DR breaks down for 
very large perturbations, when the effective perturbation 
is larger than the square root of the effective Planck's con- 
stant [T2]. However, DR remains accurate for fairly large 
perturbations, even when corresponding classical trajec- 
tories of V acc and Va PP r are completely different [7} 18} ITT]. 

The fundamental reason why quantum dynamics cal- 
culations are expensive is nonlocality of quantum me- 
chanics: Wave function ijj(t + At) at any point in space 
depends in general on ijj(t) in the whole space. There 
are many computational methods for quantum dynamics, 
but for the sake of demonstration, we consider two meth- 
ods that represent two very different general approaches. 

The first approach starts with the construction of a 
global PES, with a computational cost cn d where c is the 
cost of a single potential evaluation, d is the number of 
degrees of freedom (DOF), and n is the number of grid 
points in each DOF. Once the PES is known, dynamics 
can be performed, e.g., by the split-operator method [15] . 
In this method, the quantum evolution operator for time 
step At is approximated by 

e -iHAt/h _ e -iTAt/h e -iVAt/h + o(At 2 ), 

where H = T + V is the Hamiltonian of the system and 
T is the kinetic energy operator. Quantum dynamics 
consists of alternate kinetic and potential propagations 
(which are just multiplications in momentum and coor- 
dinate representations, respectively) and a fast Fourier 



transform (FFT) to switch the representation in between. 
The complexity of FFT is 0(N log N) where N = n d de- 
notes the dimension of the Hilbert space, so the cost of 
propagation is 0(d tn d log n). Note that the same num- 
ber of potential energy evaluations is required no matter 
how long the propagation is. This becomes an advantage 
for long-time dynamics calculations and a disadvantage 
for short-time ones. Finally, memory requirements make 
this approach feasible only for very small systems. 

In the second approach, potential energy is calculated 
"on the fly" only in the vicinity of the propagated wave 
packet. At each time step, the cost is cn d where n is 
the number of grid points in each DOF on which ip(t) is 
not negligible. Presumably, h <C n, but the exponential 
scaling with d remains. Assuming that log n < c, the 
cost of actual propagation (e.g., by FFT) is negligible to 
the cost of potential evaluation, and the overall cost of 
the dynamics is 0(ch d t). 

In the DR, potential energy and forces can also be eval- 
uated on the fly, along classical trajectories. But at each 
time step, the cost is only 0(dcn pa ths) where n pat hs is 
the number of classical trajectories used. The "hidden 
cost" in usual SC approximations is the strong depen- 
dence of ripaths on t, d, or the type of dynamics. In Ref. 
[8], it was shown rigorously that n pat h s = C(Fdr)ct -2 
where Fdr is the value of fidelity one wants to simulate, 
of dr is the error (due to finite n pa ths) that one wants 
to reach, and < C < 3. Consequently, for given Fdr 
and cr, the required number of trajectories is indepen- 
dent of t, d, or the type of dynamics! It was also shown 
that when F — > 1, which is most interesting in our ap- 
plication, C — ► 0, i.e., np a ths needed for convergence of 
Fdr becomes even smaller. The overall cost of the DR 
dynamics is 0(dcn pa ths t). In particular, there is no ex- 
ponential scaling with the number of DOF or time. 

Clearly, DR is faster than the construction of a global 
PES and than the quantum dynamics on the fly. Only 
at very long times t, in the first quantum approach, since 
the global PES is already constructed, the cost of prop- 
agation becomes dominant and one would expect that 
QM dynamics would beat DR dynamics which requires 
new potential evaluations. But even then the ratio of the 
costs of DR and QM is cn pa ,t\ is /n d log n. Assuming that 
the number of nuclear DOF is comparable to the number 
of electronic DOF, then even for the most accurate ab 
initio methods (e.g., the coupled clusters) c scales poly- 
nomially with d, and so for large enough d, DR will still 
be faster than QM dynamics. 

The DR calculation can be further accelerated by ex- 
changing roles of Vacc and Vappr- We will denote the DR 
expression defined in Eqs. by DR1 and by DR2 
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FIG. 1: Fidelity decay (top) and spectrum change (bottom) 
for D' eCO = D e ,co +0.02eV = 11.26 eV. 



FIG. 2: Fidelity decay (top) and spectrum change (bottom) 
for f?g co = Re.co + 0.02 ao = 2.15 ao. Notation as in Fig. 1. 



an analogous expression, 

/dr2 (t) := J dx Pw (x°)exp [iAS{x°, t)/h] , (4) 

AS(x°,t) = f dr AV [ql ppr (x )] . (5) 
Jo 

DR2 denotes the DR computed as an interference inte- 
gral due to action of —Ay along the trajectory of V appr . 
From definition (jTj) , it is clear that exchanging V acc and 
Vappr results in complex conjugation of /qm and has no 
effect on -Fqm- In principle, there could be an effect on 
/dr because expression Q does not have such symmetry, 
but numerical evidence presented below shows that DR1 
£3 DR2, providing further support for the approximation. 
In applications where efficiency is important, one should 
choose DR2 over DR1 since DR2 requires values and gra- 
dients of the "cheaper" PES (Va Ppr ) but only values of the 
more expensive PES (V acc ) whereas DR1 requires values 
of Vappr but both values and gradients of V^, cc ■ In applica- 
tions where calculation of DR1 is affordable, comparison 
of DR1 and DR2 results can be used as a validity test of 
the DR method since comparison with /qm will not be 
available (computation of /qm would require full quan- 
tum dynamics on V acc ). A large difference between DR1 
and DR2 results would be a sign of the breakdown of DR. 
So the requirement -Fdri ~ -Fdr2 is a necessary but not 
a sufficient condition for the validity of DR. 



One could object that since DR is an intrinsically SC 
approximation, it might suffice to estimate fidelity clas- 
sically. We explored this idea by comparing quantum 
fidelity with its classical (CL) analog, called classical fi- 
delity [14], defined in our case by 

Fch{t) := h d J rfxp t CLiacc (x)p^ Lappr (x), (6) 

where p^L ^ s the CL phase-space density evolved with the 
indicated potential to time t. For initial Gaussian wave 
packets, = p%j. Unlike DR, classical fidelity does 
not include any dynamical quantum effects, and below 
we show that indeed Fcl gives much worse results. 

To show the feasibility of our method, we have ap- 
plied it to the photodissociation dynamics of a collinear 
carbon dioxide molecule, a model that had been stud- 
ied extensively by both quantum-dynamical [THl US] and 
SC methods [T7]. Invoking the Franck-Condon principle, 
photodissociation process is described by the quantum 
dynamics of the initial state (vibrational ground state 
of the electronic ground state PES) on the dissociative 
excited PES. One can obtain the photodissociation spec- 
trum simply by taking the Fourier transform of the auto- 
correlation function C(t) :— (ip(0)\ip(t)). In future appli- 
cations, one will not be able to find V'acc(i) and /qm(^) 
due to the tremendous computational expense. Here, in 
order to demonstrate the accuracy of the method, we 
want to compare /dr(^) with /qm(^) and so we define 
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V ac c ■= Vleps to be the analytical LEPS potential for 
C0 2 [13 EH] and Kppr := V£ EPS to be the LEPS poten- 
tial with one of the parameters perturbed. Specifically, 
we change either the equilibrium bond length R Ci co or 
the bond dissociation energy -D e ,co of the CO bond. We 
imagine that we would like to obtain a spectrum corre- 
sponding to V acc but can only afford quantum dynamics 
on Vappi-. By estimating quantum fidelity amplitude /qm 
by /dr, we can determine whether we can trust the spec- 
trum computed using quantum dynamics on Va Ppr . From 
another perspective, we can use /dr to evaluate how er- 
rors in experimental values of R e ,co and D e ,co affect the 
computed spectrum. 

Figure 1 shows an example where V^ EPS has a per- 
turbed bond dissociation energy, D' e co = -D e ,co + 
0.02 eV = 11.26 eV. The figure shows that three ap- 
proaches to compute fidelity (QM, DR1, DR2) give very 
similar results. This turns out to be a small perturbation 
since fidelity remains close to unity, Fqm > 0.99. There- 
fore we should be able to trust the spectrum computed 
using ^appr(i)- This is justified in the bottom panel where 
spectra corresponding to V acc and Va Ppr prove to be al- 
most identical [20 . Unlike DR, classical fidelity (CL) 
decays much faster than its quantum analog. Judging 
by i*CL ° n lyi ° ne would conclude incorrectly that the 
spectrum computed using Va Ppr should not be trusted. 
Partially constructive interference, captured by DR, pre- 
vents quantum fidelity from decaying with the fast rate 
of classical fidelity decay. 

Figure 2 shows an example where V^ EPS has a per- 
turbed equilibrium bond length, R' e co = R e ,co + 
0.02 a — 2.15 do- The figure shows that Fdr agrees with 
Fqm and even reproduces detailed oscillations of quan- 
tum fidelity This turns out to be a large perturbation 
since fidelity falls below the value of 0.5. We therefore 
expect the spectrum computed using - 0appr(i) to differ 
significantly from the "true" spectrum computed using 
V'acc(i)- Indeed, the bottom panel shows that the spec- 
trum of Vappr is shifted and has different peak intensities 
than the spectrum of V acc . Figure 2 also shows that for 
large AV, classical fidelity starts to behave similarly to 
quantum fidelity, but is still much worse than DR. 

All calculations were performed for 600 time steps of 
0.16 fs each. Converged quantum calculations required 
n = 1024 points to discretize each CO bond length from 
to 20 a w 10.6 A, altogether using a 1024 x 1024 grid to 
represent the PES and the wave function. The DR cal- 
culation converged fully with n pat hs = 512 trajectories 
(shown in Figs. 1 and 2), but a much smaller value, 
Hpaths = 64, already gives very accurate results suffi- 
cient for our application (not shown). Both figures show 
fidelity F, rather than fidelity amplitude / which is a 
complex number containing more information. If F <C 1, 
clearly our approximate dynamics is not sufficient. How- 
ever, if F « 1, the spectrum could still be affected by a 
time-dependent phase of /. In such cases one should also 



examine fidelity amplitude. 

Our fidelity calculation for photodissociation of CO2 is 
to our knowledge the first fidelity calculation for a realis- 
tic chemical system and it is reassuring that DR remains 
valid [§]. Clearly, in chemical physics one is interested 
in systems with many more than two DOF. However, al- 
ready in the simple CO2 system, the DR calculation of 
fidelity was more than 100 times faster than the exact 
quantum calculation. We expect that in larger systems, 
much larger speedups could be achieved. 

The information required in a DR calculation is simi- 
lar to information needed in molecular dynamics (MD). 
Implementation of DR into any MD code would require a 
single addition: calculation of the action AS*. But there 
is an important difference between DR and MD: in DR, 
nuclear quantum effects are included at least approxi- 
mately, whereas in MD, even if ab initio electronic poten- 
tial is used, they are completely lost. This can be clearly 
seen in Figs. 1 and 2 where Fdr ~ Fqm 7^ -Fcl since 
Fch is basically a MD calculation of fidelity. With little 
effort, MD codes that currently compute only classical 
nuclear dynamics, could be used for evaluating accuracy 
of quantum dynamics on the same PES. 

The method was designed with the goal of determining 
the accuracy of quantum dynamics on an approximate 
PES compared to the exact dynamics. As mentioned 
above, we do not know V^ xac t and instead must use V acc . 
We predict that if one could show V acc to be much closer 
to 14 xac t than to Va Ppr , our estimate of fidelity would 
predict the accuracy compared to the exact dynamics, 
and not just compared to the dynamics on V acc . 

The DR approach is applicable to other types of per- 
turbations of the Hamiltonian than discrepancies in the 
PES. For example, DR could be used to evaluate how 
laser pulse noise affects quantum control |19j or how per- 
turbations affect quantum computation. To conclude, we 
do not claim to have found a fast way to do quantum dy- 
namics on an accurate ab initio potential. Instead we 
have found a promising method to estimate the accuracy 
of quantum dynamics on an approximate potential. 
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